ft 



REPUBLIQUE 



FRANCA1SE 




IN ST I TUT 
NATIONAL DE 
LA PROPRIETY 
INDUSTRIELLE 



REC'D 1 4 MAY 2003 



WIPO 



PCT 



BREVET D'INVENTION 



CERTIFICAT D'UTILITE - CERTIFICAT D'ADDITION 



CO 



COPIE OFFICIELLE 

Le Directeur general de Tlnstitut national de la propriete 
industrielle certlfie que le document ci-annexe est la copie 
certifiee conforme d'une demande de titre de propriete 
industrielle deposee a Tlnstitut. 



m 

I 



Fait a Paris, le 2 3 AVR. 2003 



DOCUMENT DE PRIORITY 

PRESENTE OU TRANSMIS 
CONFORMEMENT A LA 
REGLE17.1.a)OUb) 



Pour le Directeur general de i'lnstitut 
national de la propriete industrielle 
Le Chef du Departement des brevets 




Martine PLANCHE 



1NSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



SIEGE 

26 bis, rue da Saint Peterebourg 
758C0 PARIS cedex 08 
Telephone : 33 {0)1 53 04 53 04 
Taecopie : 33 (0)1 53 04 45 23 
www.bpl.fr 



» 




26 bis, rue de Safnt PStensbourg 
75800 Paris Cedex 08 

Telephone : 33 {1) 53 04 53 04 Ttlteqpie ; 33 (1) 42 94 86 54 



N° 12354*02 



\SSmm RENNES 

luEU 



I N e O'ENREGISTREMENT 
I NATIONAL ATTRIBUE PAR L'lNPI 
I DATE DE DEPfiT ATTRI8UEE 
1 PAR L'lNPI 



02 16929 
3 1 DEC. 2002 



I Vos references pour ce dossier 

\(facu!iatij) 8704 

Confirmation d'un depot parte fecooie 

Demande de brevet 
Demande de certificat d'utilite 
Demande divjslonnaire 



! 1§ 



ou demande de certificat d'utilite initiate 



BREVET D'INVENTION 

CERTIFICAT D'UTILITE 

Code de la proprias inteflectuelle - Uvre VI 

REQUITE EN DELIVRANCE 
Page 1/2 

Cetimprime esta rempllr lislblement a ^ nolrg BBM0 ^ /O1 , 

■Q mom Fct adresse l>u uemandeur ou du mandatairF 

A QUI LA CORRESPONDANCE DOIT ETRE ADReS 
Cabinet Patrice VIDON 
Le nobel - Teohnopole atalante 
2, allde antoine Becquerel 
BP 90333 

35703 RENNES CEDEX 7 



□ N° attribue par I'INPI a la telecopie 



m 
□ 
a 

N° 



Da te I i I i I | | | | 



Transformation d'une demande de TTT 
, brevet europeen Demande debrwet initiate 
S T ' TRE ^'INVENTION ,200 caries ou e SpaC e S ^ 

aitttSlSaSff* da " S « source, programme amateur 




Date 



JO DECLARATION DE PRIORITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANTERIEURE FRANCAISE 

i Nom — — '■ ' " ; 

ou denomination sociale 



Pays ou organisation 

Date l i I I If | | | 

Pays ou organisation 

Date » i I lXj j f I 

Pays ou organisation 
Date ( ) | , 1 | 



N° 




ntttntvsT 

lAfaOFtllll 
IMSU1TD1IUI 



CERTBFICAT D'UTILITE 

REQUETE Eft! DELIVRAWCE 
page 2/2 



N* D'ENREGISTREMENT 
NATIONAL ATTRIBUE PAR LINPI 



OB 540© W/ 010801 




Paiement echelonne de la redevance 

C deux tenements ) 

M REDUCTION DU TAUX 
DES REDEVANCES 



Si vous avez utilise Pimprime «Suite,> P 
mdiquez je notwbre d e pagesjptntes 
B2J SIGNATURE DU DE 
OU DU M&mBSffimE 
(Worn et $uaKte du^jgna^fre) 
P. VIDON Mang^jre (CJPi 92-1250) 



Uniquemorrtpour ies pensonnes piques effectual eWe^emes feurpropr* depot 

□ Oui 

□ Non 

Uniquement pour les personnes physiques 

□ Requ.se pour «a prendre fob pour cette invention (joMnunauis* nonunion) 
U Obtenue anterieureroent a ce de P 6t pour cette Invention Qoindre une copie de la 



VISA DE LA PREFECTURE 

mSBL. 



^^^^^^^^^ 



formulaire. 



w w 



1 



Procede et dispositif de detection de points d'interet dans une image 
numerique source, programme d'ordinateur et support de donnees 
correspondants. 

5 1. domaine de Finvention 

Le domaine de Finvention est ceiui de la detection de points d'interet, 
encore appel£s points saillants dans une image numerique. Plus pr<£cisement, 
T invention concerne une technique de detection de points d'interet mettant en 
ceuvre une approche de type ondelettes. 

10 Un point d'interSt peut etre considere comme le representant d r une region 

spatiale de F image vehiculant une partie importante de Finformation. 

Historiquement, la notion de point saillant a ete proposee dans le domaine 
de la vision par ordinateur, oh Tun des problemes majeurs consistait a detecter les 
coins des objets (d'ou le terme « saillant », utilise par la suite a titre de synonyme 

15 du terme « d'int6ret »). Plus tard, cette notion a ete elargie a d'autres 
caracteristiques des images comme les contours, les jonctions, etc. 

Dans le domaine du traitement d'images, la detection des points saillants 
correspondant aux coins des objets n'a que peu d'inttSret. En effet, les coins sont 
gen^ralement des points isoles, ne representant qu'une faible partie de 

20 rinformation contenue dans Fimage. De plus, leur detection g£nere des amas de 
points saillants dans le cas de regions textunSes ou bruit<Ses. 

DiffSrentes autres techniques ont ete proposees, s'interessant notamment 
aux points saillants correspondant aux zones de hautes frequences, c'est-a-dire 
aux contours des objets. L'invention s'applique plus pr<5cisement a ce type de 

25 technique. 

On presente ci-apres plus en detail differentes techniques connues de 
detection de points saillants. 



2 Etat actual des connaissances publiees sur le sujet (art anterieur 
le plus proche, publications, brevets) 

La detection de points saillants (aussi appeles points d'interet) dans les images est un probleme qui a suscitS 
de nombreuses recherches depuis plusieurs annees. Nous presentons dans cette section les principals 
approches retenues classiquement dans la literature. Le lecteur pourra se referer a [5] pour un etat de l'art 
plus detaill6. 

Une des premieres methodes a ete proposee par Harris et Stephens [7] pour la detection des coins. Ce type 
de points etaient alors considere comme vehiculant une quantite importante de l'information et trouvait son 
application dans le domaine de la vision par ordinateur. 

Pour definir ce detecteur, on calcule en chaque point p(x,y) de l'image /, la quantite : 

R xy = Det(M xy )-kTr(M xy ) 2 

ou M c v est une matrice definie par : 



M xy = G(a)<8> 



WW) I x (x,y)I y (x,y) 
IJx,y)I y (x,y) l;(x,y) 



ou 



<* G(a ) denote un noyau gaussien de variance a 1 ; 
*> ® denote le produit de convolution ; 

❖ l x (resp. I y ) denote la derivee premiere de / suivant la direction x (resp. y) ; 

❖ Det(M x J denote le determinant de la matrice M xy ; 

❖ Tr(M xy ) denote la trace de la matrice M v . v ; 

❖ k est une constante generalement utilisee avec une valeur de 0.04. 

Les points saillants sont alors definis par les extrema locaux positif de la quantite R xy . 

Dans [5], les auteurs proposent aussi une version plus precise du detecteur de Harris et Stephens. Cette 
version remplace le calcul des derivees de l'image / par un calcul precis des derivees du noyau gaussien. 

Le detecteur de Hams et Stephens presente ci-dessus a ete etendu au cas des imases couleurs dans [6]. Pour 
ce faire, les auteurs etendent la definition de la matrice M vv qui devient alors : 



M x y = G( a )® \ (R < +G * + % ) ( * A + O x G y + B X B V )( x .y ) 
. • l(KX y + G x G y+ B A B y )(.x,y) (R; + G;+B;)(x,y) 

ou: • . 

* d™lf^T PCC ' iVemem * Premi * reS ^ P,mS «*>*»«9» rouge, van c. bleu 

* d^f uteut n 7r PeCtiVCmem lK ^ "*» -tatofcHi- n>uge, ve« e, bleu 
Dans [10], les auteurs considered les points saillants comme les noint* i>i™*„ - . 

Au niveau k de la pyramide, le contraste au point P est defini par : 

ou G A |P; deTinit la luminance locale au point P et au niveau L- *i n /di •„ i , 

ooint P et au ni VM n h r« ^. j ? ' 1 B * P) d6finit la luminance du fond local au 

pVl^e 0„ U ^„c .^.ST^^ra^^oetr ^ ^ * ta 
fond et de'finies D ar • pyramiaes appeJees pyiarmde de luminance et pyramide du 



G k (P)= Y J ^M)G k _ 1 (M) 

M€Fih(P) 

B k (P)= XW(Q)G M (Q) 

Qe Parcntf P ) 



ou : 



Les notations Fils(P) et Parem(P) denotent les relations de hierarchie dans la pyramide oaussienne ' 

• ss^sr compte la fa - dont p ~ -~ - 

I ' 255~B t (P) J 

ZTl fixf n ° UVelle SrandCUr ' " P ° imS SaUlantS S ° nt d ' flniS Pat ,6S maXima loCaux de ^-P^rieurs a un 

^seme^Wn1ir intS P r6s t nt l initia,eraent d ans [11] est sans aucun doute le plus proche de la 

u 2l jin), et ce, a n importe quelle resolution 2'( /S-IJ A Dartir de ceu* 
coupon. u„e bie^bie de coefficient oudeteues es t co„ S<rai ,e. Ce„e bi ej/bie deLlCu^ue 




_ niveau de resolution 2'c, pour chaque coefficient ondelcttc D„f(*)* ce niveau, PemcmWeSls 
• efficients ondeieues du niveau de notation immediatement superieur 2'«Wssaires pour calculer 

Cf D,,f( n )) = {D 2 „f(k),2n<k<2n + 2p-l}osn<2 J N 

13 ^ d ' 0ndC,C " e <U - " d " <*» <~> « N denote ,a 

Ainsi, chaque coefficicm o„de,et,e Jest caicuie I partir de T> P points du signal/ Ses coefficients 

Ab , erB„/r«;; donnent la variation d'un sous-ensemble de ces T> p points. Le sous-ensemble le plus 

processus, „n coefficient D^/f selectionne a la reaction I.e. coefficient represent 2 P points du 

s^^tt- — de ,a b - ^ - z 

me^f T "° Ur ,0Ut S=/ "' i ' P ° Ur d " 9 " ° riCnlalio " 0-orizon.ale, venieale e, ob.ique) la 

comme un point saillant. a " 16renles one "'»'<™ amenent au meme pixel de I. alorc ee pixel est consider* 
Cette technique a (16 notamment utilise en indexation d'images dans [9], 

3 Inconvenients des techniques anterleures 

definition meme de point saiUant JEaoriou^n,' l« f^™ 06 i 1,a J eore da <» approches repose sur la 
s'interessaient aux coins Z 72 ^ ^".f" !. " f" erC " eUrS dans la d °™^ de la vision par ordinateur 
ete rfcemment etendu i la cotUeur danf WnZ J^F* 0 * ° < ^ eCteUr * HaiTis et Ste P he " a W <N a 
- Nation pertinente d'" 

texluiees, ces points vont etre eparpilles dans l'eZSTJ i ' CaS d lmages fakement 

l image. Dans le cas d'images S urles ou E w donneront pas une representation satisfaisante de 
'--e.donneron^^ °"°»' <™ ™* da " a 

d a nttrr, h ?u^mlTceX^ P°« le traitemen, 

texlurees ou bmitees. PP f ™ d " meme d<5faul ,a P^adenle dans le cas de regions 

™"i^^:T*«JT11r " N ' Sebe »" «" <a plus 

primord,ale d' une taSJ^SNS e^ fSS^TT *" *" C0I "° Ure "I*"™ rfafcnJL 

app n e « ro t *- ^^P^s^^^rf"" visuel homain - caire 

' pa^~^c^ """T^ ProP ° Si ^ Lo " ta « tb ° ^ « ~e 
«3,h ^coefficjems pour chaque niveau de resolution 2>et pour une image carree Notre 

~? P r iCxicr^" rr de ,a base <,, ° nde,ei ' e: * - J: ^ ^ : 

d'ondelene avec un support ™TL * , a n ^ n ° US PCn " e ' de pOUVoir u,iliaer * a bases 
utilisantled.teoteurdeKrsrutiltnTi^ 



» 



""ZD' 



La methode de Loupias ct Sebe considered lcs sous-bandes independamment les unes des ,„•«.« ™ 
qui lcs aminent a dctccter en priorite les points de gradient maximal dans toutcs les directLt i e 
les corns). De noire c6td nous fusionnons rinfonnation contenue dans les differentes souTbande ce 
qui nous permet de ne privilegier aucune direction particuliere. 

4 Le but de I'invention et !e resuitat technique obtenu par 
['invention 

Le but de ^invention concerne la detection des points saillants d'une image /. Ces points correspondent aux 
pixels de / appai enant 4 des regions de haute frequence. Pour ce faire, nous proposons de^^aser sur la 

cen°e"/orie° ^ [1K2][33 - N ° US rem ' 0y0nS k IeCl£Ur " rAnneXe A P~ brL p" S enSr d e 

La transforms en ondelettes est une representation multi-resolution de 1'image qui permet d'exprimer 

1'image aux differentes resolutions 1,1, etc.. Ainsi, a chaque niveau de resolution 2'(j&-l). la- 

transformed en ondelettes represente 1'image /, de taille nxm = 2*x2' (k,l& Z)\ sous la forme : 

❖ d'une image grossiere \,I ; 

❖ d'une image Z>',/de details representant les hautes frequences verticales (i.e. les contours 
horizontaux) ; 

❖ d'une image Dj/de details representant les hautes frequences horizontales (i.e. les . .contours 
veilicaux) ; 

❖ d'une image Z)J/de details representant les hautes frequences diagonals (i.e. les coins) ' 

clZZl f ^ taillC ^ X2 ' +/ • La RgUre 1 ilIUStre ce ^ e de -presentation. . I 

Chacune de ces trots images est obtenue a parlir A^I par un filtrage suivi d'un sous-echantillonnage d'un 

facteur deux dans chaque direction comme le montre la Figure 2. U est a noter que Ton a A / = / V- 

L mvention consiste a choisir, en premier lieu, une base d'ondelettes et un niveau de rlblution^nimal 

r (> <-D. Une fois la transformation ondelette effectuee, nous proposons de parcourir chacune des trois 

usages de details Dy I . D;J et Dj/afin de construi re une arborescence de coefficients ondelettes Cette 

P^r d Tl e l b f 6 7 raPPr ° Che Zer0tree C4L iniUa,6nient Pr0p ° s6e P ° Ur Ie Coda ^ d Elle nous 

coXi : f r^ir ?<:tti? ll : 2k ~f- refl6tant nm ~ de ^ 

id resojuuon z { i S-l). Ainsi, un coefficient ayant une saillance imnnrtiTifp 

S" 1 \"" *f °" de ' PrtSentan ' d « ha " KS "*"«• E " - ^ncta, o„«e«e d"~ 
a la resoiuuon 27r<-l, con-espond a u» contour de .'image ^./suivan, une direction 
pamcuhere (honzontale, verlicale ou oblique). L'approche Zerotree nous indique que chacun des 
SSf deT 4 „' a r& ° MOn 2 Ws P0 " d 4 — de Wile 2~ X2~dans rtaage , A 

-'X2 riTlf I n T 1 ™"- '' inVeMi0n Pr ° P0Se ° ne m ^ 0de P ermeMa " 1 * choisir paSni let 
- xz pixels de/.le pixel le plus representatif de cette zone. 

ra^^ P ° temielleS ' 13 d ' leCli0n dC P ° imS SailIantS danS les ™^ P« •» de 

❖ Pour le tatouage d'images. Dans ce cas, les points saillants donnent des indication auanf h T* 
9 localisation possible de la marque afin de garantir sa robustesse ■ indication quant a la 

V sZ!l" d ?~ a v 0n En d ' tCCtant Un n ° mbre flXe d ^ P oi " ts ^illants, on pent en deduire une 

. ====i^^ — - ' — > - P- 

Fimaoe ^Td't,? 88 " 8 - ^ P ° inlS Sail,antS ^espondant aux hautes frequences de 
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L'invention a done notamment pour objectif de pallier les differents 
inconvenients de l'etat de Tart. 

Plus precisement, un objectif de l'invention est de foumir une technique de 
detection de points saillants correspondant a une haute frequence, et ne 
privilegiant aucune direction particuliere dans 1 'image. 

Un autre objectif de l'invention est de foumir une telle technique, qui 
necessite un nombre reduit d'operations, par rapport aux techniques connues. ' 

Notamment, un objectif de l'invention est de foumir une telle technique 
permettant d'utiliser des bases d ondelettes avec un support de taille importante. 

Ces objectifs, ainsi que d'autres qui apparaitront plus clairement par la 
suite, sont atteints a l'aide d'un precede de detection de points d'interet dans une 
image numerique source, ledit precede mettant en ceuvre une transformation en 
ondelettes associant a une image source des coefficients ondelettes, pour au moins 
deux niveaux de resolution : 

- une image grossiere ; 

- une image de details representant Jes hautes frequences verticals ; 

- une image de details representant les hautes frequences horizontals ; 

- une image de details representant les hautes frequences diagonales, 

un point d mt6r€t etant un point associe a une region de 1'image presentant une 
haute frequence. 

Selon l'invention, ce precede comprend les etapes suivantes : 

- application de ladile transformation en ondelettes a ladite image 
source ; 

- construction d'une arhorescence unique, a partir des coefficients 
ondelettes de chacune desdites images de details ; 

- selection d'au moins un point d'interet, par analyse de ladite 
arborescence. 

De facon avantageuse, ladite etape de constmction d'une arborescence 
repose sur une approche de type Zcrnlree. 

Ainsi, preferentiellement, chaque point de 1 -image grossiere de resolution 



minimum et la racine d'un arbre a laquelle on associe trois nceuds fiis 
correspondant £ chacun des coefficients ondelettes des trois images de details 
localises au mSme emplacement, puis on associe recursivement, a chaque noeud 
fils d'un niveau de resolution donne, quatre noeuds fils formes par les coefficients 
5 ondelettes de Timage de details de meme type et du niveau de resolution 
precedent, et associee a la region correspondante de I'image source. 

Selon un aspect avantageux de Finvention, ladite etape de selection met en 
oeuvre une etape de construction d'au moins une carte de saillance, affectant 
auxdits coefficients ondelettes une valeur de saillance representative de son 
10 interet Prefdrentiellement, on construit une carte de saillance pour chacun desdits 
niveaux de resolution. 

De fa9on avantageuse, pour chacune desdites cartes de saillance, on 
fusionne pour chaque valeur de saillance les informations associees aux trois 
coefficients ondelettes correspondant aux trois images de detail, de fa?on a ne 
15 privilegier aucune direction dans Timage. 

Selon un aspect preferentiel de Finvention, une valeur de saillance d ? un 
coefficient ondelette donne d'un niveau de resolution donne prend en compte la 
ou les valeurs de saillance des coefficients ondelettes descendant dans ladite 
arborcscence dudit coefficient ondelette donne. 
20 Pr6f£rentiellement, une valeur de saillance est une relation lin£aire des 

coefficients ondelettes associes. 

Dans un mode de realisation particulier de Finvention, la valeur de 
saillance d'un coefficient ondelette donne est calcuiee & partir des equations (I) 
precisees par la suite. 

25 Dans ces equations, le parametre a k peut par exemple valoir -1/r pour 

toutes les valeurs de k. 

Selon un autre aspect preferentiel de r invention, ladite etape de selection 
comprend une etape de construction d'une arborescence desdites valeurs de 
saillance, repose avantageusement sur une approche de type Zerotree. 

30 Dans ce cas, ladite etape de selection comprend avantageusement les 



a 



e tapes de : 

- tri par ordre decroissant des valeurs de saillance de la carte de saillance 
correspondant a la resolution minimum ; 

- selection de la branche presentant la valeur de saillance la plus elevee pour 
5 chacim des arbres ainsi tri£s. 

Selon un aspect preTerentiel de 1 'invention, ladite dtape de selection de la 
branche presentant la valeur de saillance la plus elevee met en ceuvre un parcours 
de l'arbre correspondant a partir de sa racine, et une selection a chaque niveau de 
l'arbre du noeud fils pre'sentant la valeur de saillance la plus elevee. 
10 Comnie deja mentionne, l'invention permet d'utiliser de nombreuses 

transformations en ondelettes. Dans un mode de realisation particulier, on met en 
oeuvre la base de Haar. 

Dans un mode de realisation particulier, on choisit comme niveau minimal 
de resolution 2 4 . 

15 Le proceed de l'invention peut par ailleurs comprendre une etape de calcul 

d'une signature d'image, a partir d un nombre predetermine de points d'interet de 
ladite image. 

Ladite signature peut ainsi notamment etre utilisee pour 1'indexation 
d'images par Ieur contenu. 

20 Plus generalement l'invention trouve des applications dans de nombreux 

domaines, et par exemple pour : 

- le tatouage d'images ; 

- 1'indexation d'images ; 

- la detection de visage(s) dans une image. 

25 L'invention concerne egalement les dispositifs de detection de points 

d'interet dans une image numerique source mettant en ceuvre le precede tel que 
decrit ci-dessus. 

L'invention concerne encore les programmes d'ordinateur comprenant des 
instructions de code de programme pour l'execution des €tapes du procede de 
30 detection de points d'interet decrit ci-dessus, et les supports de donnees 
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humeriques utilisables parun ordinateurportantun tel programme. 

D'autres caractenstiques et avantages de 1'invention apparattront a la 
lecture de la description suivante d'un mode de realisation preTerentiel, donne a 
titre de simple exemple illustratif et non limkatif, et des dessins annexes parmi 
lesquels : 

la figure 1 illustre le principe de l'analyse multi-resolution dune 
image I par transformation en ondelettes ; 

la figure 2 presente de facon schematique une transformation 
ondelette ; 

la figure 3 est une representation d'une arborescence de 
coefficients ondelettes selon 1'invention ; 

la figure 4 presente un exemple de cartes de saillance, et les arbres 
de saillance correspondants ; 

la figure 5 illustre la saillance d'une branche de l'arbre de la figure 
4; 

, les figures 6a et 6b illustrent des resultats experimentaux du 
procede" de Tinvention, la figure 6a presentant deux images 
originales et la figure 6b les points saillants conespondants ; 
la figure 7 illustre un procede d r indexation d'images mettant en 
ccuvre le procede de detection de 1'invention. 




5 Identification des elements techniques essentiels de I'invention 

5. 1 Tranformation ondelette 

echanUJIonnea256vaIeurs.DepIus,noussupposonsque H = 2*(*eZJ etque m = 9'// e ZJ 
Si nous notons / 1'image originate, nous avons alors : 

I [0,m]x[cu]->[0,255] 
' (x.y)^> I(x,y) 

icTTZlT^ * " r U ° n , * ' a "T fomta ° nde,e " e * 7 1— une -Potation mute-resotaion 
de / A chaque mveau de r=so lull o„ 2'OS-U. la represent*™ de /« donnee par une image grossiere 

A,,/ et par trois images de details D [ I D 2 let r> 3 / nh*™~- ^ * , . 

s J5> » ^ 2 / y et - Chacune de ces images est de taille x Ce 

processus est illusire sur la Figure 2> " 

La transformation ondelette necessi.e le choix d'une fonction echelle <* x) ainsi que le choix d'une 
fonc ion ondelette V(x) A partir de ces deux fonetions, on derive un filtre echelle //et un filtre onTelette C 
dont Ies reponses impulsionnelles respectives h et g sont definies par : onaejene G 

M n ) = ((# JH (u),tf u - « ;) Vn e z 

S( » J = (^ 2 -. f u ), <p(u-n ))\fn e Z. 

Notons respecti vement H et G Ies filtres miroirs de // et G (i.e. />„ ; = h( - n ; et I, „ , = _„ . > 
On peut alors montrer [IJ (cf. Figure 2) que: 6 

❖ ,4/ peut etre calculee en convoluant ^,.,/avec tfdans Ies deux dimensions et en sous- 

echantilionnant d'un facteur deux dans les deux dimensions ; 
*> D\, I peut etre calculee en : 

1. convoluant A^I avec H suivant la direction y et en sous-echantillonnant d'un facteur deux 
suivant cette merne direction ; 

" H wlT 1 h r6mltat de r " tape 1} 3VeC G Suivant ,a direction * et en sous-echantillonnant 
dun facteur deux suivant cette raeme direction. wnamuionnant 

*> D;, I peut etre calculee en : 

1. convoluant A^l avec G suivant la direction y et en sous-echantillonnant d'un facteur deux 
suivant cette raeme direction ; 

2 ' dZtT H re ' SUUat dC r ' lape 1} aV£C '^ Suivant ,a direction', et en sous-echantillonn.m 
dun facteur deux suivant cette meme direction. ctnanujionnani 

* D 3 2 , J peut etre calculee en : 
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1. convoluanl ^/ avec Gsuivanl Ja direction y et en sous-echanlillonnant d'un facteur deux 
siiivant cette meme direction ; 

2. convoluant lc resullat de Fetape 1) avec Csuivant la direction X et en sous-echantillonnant 
d un facteur deux suivant cette meme direction. w-namijionnnnt 

5.2 Construction de I'arborescence des coefficients ondelettes 

Une Ibis la transformation ondelette effectuee jusqu'a la resolution T( r <,-!), on dispose : 

❖ D 'une image approximee A,,/ ; 

❖ De trois images de details D\, / , Z>> / , D\, I par niveau de resolution 2> avec j=-l r 

On constant alors une arborescence des coefficients ondelettes en se basant sur !a technique Zerotree T41 
Les arbres sont consu-uils de la facon suivante (cf. Figure 3) : recnnique Zerotree [4]. 

❖ Chaque pixel p(x,y) de l'image A,, / est la raclne d'un arbre ; 

❖ Chaque racine p(x.y) se voit attribuer trois ncuds fils designes par les coefficients ondelettes des trois 
images de deta.ls D*. I (s=l,2,3) localises au meme emplacement (x,y) • 

❖ Du fait du sous-^chantillonnage d'un facteur deux effectue par la transformation ondelette a chaque 
changement de resolution, chaque coefficient ondelette c$(x.y) (s=l,2,3) correspond a une zone de 
taille 2x2 pixels dans l'image detail correspondante a la resolution 2^. Cette zone est localise en 
( — v,2_y ) et tous les coefficients ondelettes y appartenant deviennent les nceuds fils de a* (xy) 

De facon recursive, on constn.it I'arborescence dans laquelle chaque coefficient ondelette of (x y)(s=l 2 3 
ct 0>u>r) possede quatre no,uds fils designes paries coefficients ondelettes de l'image D>J localises dans 
ia region situee en ( 2x,2y ) et de laille 2x2 pixels. 

Une fois I'arborescence construite, chaque coefficient ondelette ^(x,y K s=l,2, 3) correspond a une region 
de taille 2~ r x2 -r pixels dans l'image detail £>*.,/ . 

5.3 Construction des cartes de saiilance 

t«T* ^ I ; i arb ° re ! CenCe ° blenue P ar ] ' 6ta P e P^dente, nous proposons de construire en ensemble de -r 
canes de saiilance (,.e. une carte de saiilance par niveau de resolution). Chaque carte de saiilance S 0=- 

un colfficSnt 1 l^T° e d6S COefficiems ondeleltes P^ents a la resolution correspondante 2' . Ainsi] plus 
Sl^^SS^" JUSe imP ° rtam aU SenS * ! ' in "" *™ vehicule, plus sa valeur^ 

suivant l^nlfH^^T ^f^*? 1 ° nde,ette privi, ^ ie une direc «on (horizontals verticale ou oblique) 
aTrec^on Z S £ * * J Cependant, nous avons choisi de ne privildoier ^aucune 

ondetnes « , t V"", T" *"* rinfo ™<- contenoe dans les 'trois ££££ 

ondelettes ^fx.,^fx.>A<f Jquelque soit le niveau de resolution V et quelque soit la 
localisation (x.y) avec 0 < x < 2 k * J etO<y< 2 l *> 
Chaque carte de saiilance S„ est de 'taille 2 k¥} X2'*' 

resolution 2 est donnee par la relation recursive suivante : • 
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Equation 1: expression de la saillance d'un coefficient 



♦ (s=, ,3) dene* „ va.eur m a, iraol e d=s coefficient o„de Ict ,e s dans ,- taage dftail 

♦ "Jo**, <l,)permet de regler Importance des coefficients de saiUance suivant ]c niveau de 
resolution. Best a noter que Ton a =7. 

♦ II est important de no.er que les valeurs de saillance sont normalises i.e. 0 < S ( x y ) < 1 

2ee « -at i0 ; d ; coeffjcienls 

comme des pixels de y im&t appJtenan^ * "de , COnsid6 ' ons les Points saillants 

</a-,vK S =1,2 3) eleve a Ta reso^Z g * fr£qUenCe - ° r « U " coe «^"t ondelette 

2 ,< o M w; etexe a Ja icsoJut.on 2'denote une zone de haute frequence dans l'im«i» e >4 i, 

resolution 2' et corrcspondant a Porientation du contour 

Aim,, la formulation de la saillance d'un coefficient donnee dans liquation 1 se justifie. 

5.4 Choix des points saillants 

conanme 2'"- arbres de coefficl T* v ^T™" COefficieMS <»*"«'«. ™»* pouvons 

Afin de localiser les points les plus saillants dans 7, nous procedons : 

1. auntri par ordre decroissant des 2— valeurs de saillance present dans S • 

2. a la selecuon de la branche de saillance maximale de chacun des 2— arbres ainsi tries 

imponante (cf. Figure 5). Nous obtenSlf une^st^S d^sSnce- ^ * '* P, " S 

f . Wt Jfox{f r ^2x w + „-,2_y,_, + v ;,0<„ <l,0<v<l} 

A partir des branches les plus saillantes de chaque arbre le oixel H, / ,k ■ - 

repre-sentaiif de la branche est localise" en Ox o„ , P J de 1 cho,sl comme elam ^ plus 
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Dans la pratique, on parcourt uniquement un sous-ensemble des 2 k *>^ arbres. En effet, pour de nombrcu^ 
apphca uons, on recherche un nombre fixe „ de points saillants. Dans ce cas, il convient de n p^oS aue 
les /» arbres ayant les racines les plus saillantes. paicounr que 

6 Description detaifiee d'au moins un mode particulier de 
realisation 

Dans cette section, nous utilisons les elements techniques presented dans la section precedente dont nous 
fixons les parametres necessaires afjn de decrire un mode particulier de realisation. 

6. 1 Choix de fa transformation ondelette 

Comme nous I'avons mentionne en section 5.1, nous devons en premier lieu choisir une base d'ondelettes et 
le m veau mimmal de resolution 2 r (r<~l). 1 

Pour ce mode particulier de realisation, nous proposons d'utiliser la base de Haar et r=-4 
La base de Haar est definie par : 



(lsiO<x<l 
0 sivon 



pour la fonction echelle, et par 



lsiO<x<- 
2 

- 1 si — < x < 1 



0 sw on 
pour la fonction ondelette. 



6.2 Construction de i'arborescence des coefficients ondeiettes 

section^ 6taPe ' aUCUn Param ' tre "'^ reqUiS - U prcCeSSUS eSt donc conf ^e a ce qui est decrit dans la 

5.3 Construction des cartes de sa ilia nee 

IZ^u^r^ Ch0iSk kS a ^ k ^ de regler Importance 

accoide aux coefficients de saillance suivant le niveau de resolution auquels ils appartiennent. 

Dans ce cas particulier de realisation, nous proposons d'uuliser a k = — Vk e [r,-l]. 

6.4 Choix des points sailiants 

Cette elape ne requiert aucun parametre. Le processus est done conform a ce qui est decrit dans la section 

6.5 Resuitats experimentaux 
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6.6 Exemple d'application 

Parmi les applications potentielles listees dans )a section 4, nous presentons dans cette ser.inn iwr 

des points saillants pour 1'indexation d'images fixes par lecontenu " 1 Ut,hSatl0n 

6.6. 1 But de indexation d'images 
L'indexation d'images par le contenu permet de retrouver, parmi une base de donnees d'imacres Iln 
ensemble damages visuellement similaires a une image donnee appelee image reaue e Pour ce fafre d^ 

To* to d,fficu lt e de r,„dexa,io„ d'images consist aic a deten^er del descriptors e« des distances 
S.6.2 Descripteurs base's sur les points saillants d'une image 

D(R,I j ) = Z W ,S j (f i ),j = l,...,N 
ou denote le nombre d'images dans la base de donnees et S/fi) est definie par : 

penmen, d e mo du,e r Fi m pLa„ce ' dea lescrfp.eu^ ,2 U „s ^'at^ 
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Annexe A : Survo! de la theorie des ondelettes 
A.1 Introduction 

«^™%T£T^ ] ^^:™ r fon r, (courbe - «* » «*— 

^-^dadaailsp^ 

detaUs dune function a ,a Solution ,,.som definis conmre ,a differenee d'info.ruation entTson 
approximate a la resolution r,. et son approximation a la resolution r M . 
A.2 Notations 

^^J^^^^^- ~soIuUo n , nous p^ons dans cette 

Z ^T,!r b,eS , deS n ° mbreS emierS et r6eJs sont ^Peclivement not* Zet 7?. 

. L r^denote 1 espace vectorieJ des fonctions unidimensionneJIes^J mesurables et intfgrabJes 

(/c -v a *r -*;) = jT/r « ;^r« . 

* Pour f(x)eI?(R)« ^^^V^;jaconvo]uuonde/ W et 5WeS tdefiniepar: 

/*Sf*J = $~J(u)g(x-u)du. 

❖ Zrf ^ ^denote J'espace vectorieJ des fonctio"ns ^ de deux variables mesurables et integrates 
- Pour f (x , ,;e ,et A le produit scalaire de/,,,, et ^; est de*ni p" 

(/f*.;y;.tf *.yj)= JT £f(u,v) g (u, v )diid v . 
A.3 Proprietes de I'analyse multi-resolution 

Les proprie^s attendues de A,, sont les suivantes : 

1. A v est un operateur lineaire. Si ,,repre-sente Fapproximation dcff.i) a la resolution 2' alors 

A f( X )ne doit pas ft* modifie lorsqu'on 1'approxime de nouveau a la resolution 2'. Ce principe 
s ecnt ^ . ^ . ^ et montre que rop , raleur K> ^ ^ op&ateur ^ d ^ ^ 

vectonel I/, cZrf*;. Cet espace vectoriel pent etre interpret comme 1'ensemble de toutes les 
approbations possibles a la resolution V des fonctions de L 2 ( R) 

2. Routes les approximations possibles de/ W a ,a resolution V, ^,est la plus sinulaire a 
/fa;. L operateur A 2 , est donc une projection orthogonale sur V,, . 

3. L'approximation d'une fonction a la resolution ^ +, cont,Vnt i- # 

caJculer la meme fonction a la r&olution hO^Lo'lZT ' 1 f 0 ™ 1 ™ n<kessaire P our 
suivante : resolution mteneure 2 . Cette propnete" de causalite induit la relaUon 

Y/eZ.V^cl/,,,. 




4. L'operation d' approximation est la meme a toutes les nfenint.™* t „ 

EST "* - - - -= - - &ffi££,2SK 

Vj£Z,f(x)e V,, 2.v J e V^, . 

5. Lorsque Ton calcule une approximation de fix) a la resolution 2' une oartie de 1'infn <• 
contenue dans fix) est perdue. Cependant, lorsque la resolution te'nd ve^Mnfini 1 f° f nTiahon 
approximee doit converger vers la fonction fix) oriainale De ! ii! f , \ 3 f ° nCtl0n 
tend vers zero, la fonction approve cJa££S^^ 

Tout espace vectoriel f V„ ; /s2 qui satisfait IWble de ces proprietes est appele ^ 
resolution de L~(R). 

A.4 Analyse multi-resolution d'une fonction unidimensionnelle 
A.4. 1 Recherche d'une base de v 

Nous avons vu dans la section A3 que Topfateur ^approximation A, ea« une projection otlhogonaie sur 
1 espace vec.one! V„, Afin de caracteriser numeriquemen, cet operateur. nous devons trouver une base 
orthonomale de V,, . 

V^etant un espace vectoriel contenant les approximations de fonctions de l}( R) a la resolution 2' toute 

:r ypeut 6ue vue comme un vecteur * 2> — ■ * ~ f - *« — * 

^,t S L^f meS r ri T PaUX ^ 13 th '° rie ° nddetteS St * U,e *™ «** une fonction unique 
^ W ' 3PPelee f ° nCt,0n 6chdle ' 4 P artir de Ia <^ on peut definir 2> fonctions de base l^Z 
V 2J par dilatation et translation de<Z>(x): ■ ' ' ' 

*/r = or 2'*-;;, = o,-.,2' -1 

sur ies 2 foncnona de base *>(*>■ Cette operation eonsiste a calcuier .e produi, scalaite de f M avec 
chacunedes 2' fonctions debase */("x; : JW 

= lL\f( u )M2 j u-k))b(2 j u-k). 
: 1 ^ ^ PEUt *" rtdUit * COnVOh " i0 " "« h fil - ~* 

Comme <b( x ) est un filtre passe-bas A f r>eut Stn» int^^i-s ^ 

eehantiiionnageuniforme. ' «re mterprete comme un filtrage passe-baa suivi d'ut, soua- 

A.4.2 Construction de /'analyse multi-resolution • 

Soh A fv» n ■ ■ , f ^ ; echantllJons - La reso ^tion maximale devest alors n. 
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T ,a proj T des 2,fonc,ions de base * 1,x>de v * ™ v ~ • - — *» 

peut etre obtenue en convoluant /I /" aver \e fiitr* r,o<.<. a u„ _* 
" . . 2'*' 7 ec le mtre Passe-bas correspondant a Ja fonction echell^ Pt * n 

sous-e-chantillonnant le resultat d'un facteur 2 : toncuon echelle et en 

^ 2 y/C" ; = 2> A- - 2u )A ,f(k ),0<u<2 J -l 

k=0 

avec ^ n J = (ty 2« ®(u-n )), Vn g Z . 
La fonction detail 

Comme mentionne dans la propriete C5) de la section A i' A „^„ t - 

fix) a une resolution 9' a L*; aelasectl0n A - 3 > 1 operation qui consiste a approximer une fonction 

approximee A,,*,/ COnnaJSSance de 4,/permet de reconstruire parfaitement la fonction 

c es, d,re 2 foncons de base. Un autre theoreme important de <a meorie des onde.eues stipule qu'a 
pamr d une function eeheiie « x) . „ est poss.bie de definir * functions de base de ^ . Ces fonc[j J ^ 
base*, W son, obtenues par diiatadon et tfanslaUon d'une fonction appe.ee fonction ondelette • 

*/fx;=¥f2'i-/Ai=a...^/_ 1 . 

LI" 14 "" fa9 °" P ° Ur 13 C °" Sln,Cti0n * ''WroximaMon „,,/ , on pent rnontrer que D f pent 6,re 

D„f=( f(u)'f<-2>u)Xk),keZ. 
A.4.5 Extension a I'analyse multi-resotutlon rf e fonction bi-dimensionneUes 

mani4rc ^ mui — - — ~ 

Pour ce faire, on utiiise to, meme theoremes que ceux utilises precedemmen, Ainsi, si .'on noic V respace 
vectoriel des approximations de L'(K> )i !a resoiu.ion 2'. „ peu, montrer que Ton pen, louver une b,se 
o n honorTnaledet/ j ,e„dilatante,transIa.a„tunefonctionfchelie*fx.v; ei =,^-; : 

*/r x.y ) = <D(2lx-i,2>y-j>.( i.j >sZ\ 
Dans le cas particulier des approximations separables de L 2 ( Jt 1 ) on a W«J=*fr«J) - ^ . 

montrer que si Ton note 47 v II, ft, ,• \, "" ' " ' $ b,d,ma »"»<^l °" pen, 

fonctions denniespar "°" "" CM 4 '° f ° nCli °" &hclte *» .es „ois 



REVINDICATIONS 

1. Procede de detection de points d'interet dans une image numerique source, 
ledit procede mettant en ceuvre une trans formation en ondelettes associant a une 
image source des coefficients ondelettes, pour au moin.s deux niveaux de 

5 resolution : 

- une image grossiere ; 

- une image de details representant les hautes frequences verticales ; 

- une image de details representant les hautes frequences horizontales ; 

- une image de details representant les hautes frequences diagonales, 

10 un point d'interet etant un point associe a une region de Timage presentant une 
haute frequence, 

caract£rise en ce qu'il comprend les etapes suivantes : 

- application de ladite transformation en ondelettes k laditeVimage 
source ; 

15 - construction d'une arborescence unique, a partir des coefficients 

ondelettes de chacune desdites images de details ; 

- selection d'au moins un point d'interet, par analyse de ladite 
arborescence. 

2. Procede selon la revendication 1, caracterise en ce que ladite etape de 
20 construction d'une arborescence repose; sur une approche de type Zerotree. 

3. Procede selon la revendication 2, caracterise en ce que chaque point de 
r image grossiere de resolution minimum et la racine d ? un arbre h laquelle on 
associe trois nceuds fils correspondant k chacun des coefficients ondelettes des 
trois images de details localises au meme emplacement , 

25 puis en ce qu'on associe r£cursivement, h chaque nceud fils d ? un niveau de 
resolution donne, quatre nceuds fils formes par les coefficients ondelettes de 
Timage de details de meme type et du niveau de resolution precedent, et associee 
& la region correspondante de Timage source. 

4. Procede selon Tune quelconque des re vesications 1 & 3, caracterise en ce 
30 que ladite 6tape de selection met en ceuvre une <£tape de construction d'au moins 



une carte de saillance, affectanl auxdits coefficients ondelettes une valeur de 
saillance representative de son interet. 

5. _ Procede selon la revendication 4, caracterise en ce que Ton construit une 
carte de saillance pour chacun desdits niveaux de resolution. 

6. Proced6 selon Tune quelconque des revcndications 4 et 5, caracterise en ce 
que, pour chacune desdites cartes de saillance, on fusionne pour chaque valeur de 
saillance les informations associees aux trois coefficients ondelettes correspondant 
aux trois images de detail, de facon a ne priviiegier aucune direction dans l'image. 

7. Procede selon l'une quelconque des revcndications 4 a 6, caracterise en ce 
qu'une valeur de saillance d'un coefficient ondelette donne d'un niveau de 
resolution donne prend en compte la ou les valeurs de saillance des coefficients 
ondelettes descendant dans ladite arborescence dudit coefficient ondelette donne. 

8. Procede selon l'une quelconque des revendicalions 4 a 7, caracterise en ce 
qirune valeur de saillance est une relation Iineaire des coefficients ondelettes 
associes. 

9. Procede selon la revendication 8, caracterise en ce que la valeur de 
saillance dun coefficient ondelette donne est calcuiee a partir des equations (1). 

10. Procede selon la revendication 9, caracterise en ce que le parametre a k 
vaut -1/r pour toutes les valeurs de k. 

11. Procede selon l'une quelconque des revendications 4 a 10, caracterise en ce 
que ladite etape de selection comprend une etape de construction d^iine 
arborescence desdites valeurs de saillance. 

12. Procede selon la revendication 11, caracterise en ce que ladite etape de 
construction d'une arborescence desdites valeurs de saillance repose sur une 
approche de type Zerotree. 

13. Procede selon l'une quelconque des revendications 1 1 et 12, caracterise en 
ce que ladite etape de selection comprend les etapes de : 

- tri par ordre decroissant des valeurs de saillance de la carte de saillance 
correspondant a la resolution minimum ; 

- selection de la branche presentant la valeur de saillance la plus eievee pour 



chacun des arbres ainsi tries. 

14. Procede selon la revendication 13, caracterise en ce que ladite etape de 
selection de la branche pr6sentant la valeur de saillance la plus elevee met en 
oeuvre un parcours de Farbre correspondant h partir de sa racine, et une selection a 

5 chaque niveau de Farbre du noeud fils presentant la valeur de saillance la plus 
elevee. 

15. Procede selon Tune quelconque des revendications 1 a 14, caracterise en ce 
que ladite transformation en ondelettes met en oeuvre la base de Haar. 

16. Procede selon Fune quelconque des revendications 1 a 15, caracterise en ce 
10 que le niveau minimal de resolution est 2' 4 . 

17. Procede selon Fune quelconque des revendications 1 k 16, caracterise en ce 
qu ? il comprend une etape de calcul d'une signature d'irnage, a partir d % un nombre 
predetermine de points d ? interet de ladite image. 

18. Procede selon la revendication 17, caracterise en ce que ladite signature est 
15 utilisee pour Findexation damages par leur contenu. 

19. Application du procede de detection de points dlnteret dans une 'image 
numerique source selon Tune quelconque des revendications 1 a 18 & au mbins un 
des domaines appartenant au groupe comprenant : 

- le tatouage d'images ; 
20 - Findexation d'images ; 

- la detection de visage(s) dans une image. 

20. Dispositif de detection de points d'interet dans une image numerique 
source, mettant en oeuvre une transformation en ondelettes associant a une image 
source des coefficients ondelettes, pour au moins deux niveaux de resolution : 

25 - une image grossiere ; 

- une image de details representant les hautes frequences verticales ; 

- une image de details representant les hautes frequences horizontals ; 

- une image de details representant les hautes frequences diagonales, 

un point d'interet etant un point associe a une region de Fimage presentant une 
30 haute frequence, 



caracterise en ce qu'il comprend : 

- des moyens d'application de ladite transformation en ondelettes a 
ladite image source ; 

- des moyens de construction d'une arborescence unique, a partir des 
coefficients ondelettes de chacune desdites images de details ; 

- des moyens de selection d'au moins un point d'interet, par analyse de 
ladite arborescence. 

21. Programme d'ordinateur comprenant des instructions de code de 
programme pour 1 'execution des etapes du procede de detection de points 
d'interSt dans une image numenque source selon l'une quelconque des 
revendications 1 a 19. 

22. Produit programme d'ordinateur comprenant des instructions de code de 
programme enregistre" sur un support utilisable dans un ordinateur, 
comprenant des moyens de programmation lisible par ordinateur pour la mise en 
oeuvre d'une transformation en ondelettes associant a une image source des 
coefficients ondelettes, pour au moins deux niveaux de resolution : 

- une image grossiere ; 

- une image de details representant les hautes frequences verticales ; 

- une image de details representant les hautes frequences horizontales ; 

- une image de details representant les hautes frequences diagonales, 

un point dmter€t etant un point associe a une region de l'image pr6sentant une 
haute frequence, 

caracterise en ce qu'il comprend : 

- des moyens de programmation lisible par ordinateur pour effectuer 
f application de ladite transformation en ondelettes a ladite image 
source ; 

- des moyens de programmation lisible par ordinateur pour effectuer la 
construction d'une arborescence unique, a partir des coefficients 
ondelettes de chacune desdites images de details ; 

- des moyens de programmation lisible par ordinateur pour effectuer la 



selection <Tau moins un point d'interet, par analyse de ladite 
arborescence. 

23. Support de donnees numeriques utilisable par un ordinateur, caract£rise en 
ce qu'il comprend des instructions de code de programme d'un programme 
d'ordinateur selon Tune quelconque des revendications 21 et 22. 
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